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Opacity of electromagnetically induced transparency for quantum fluctuations 
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We analyze the propagation of a pair of quantized fields inside a medium of three-level atoms in 
A configuration. We calculate the stationary quadrature noise spectrum of the field after propa- 
gating through the medium, in the case where the probe field is in a squeezed state and the atoms 
show electromagnetically induced transparency (EIT). We find an oscillatory transfer of the initial 
quantum properties between the probe and pump fields which is most strongly pronounced when 
both fields have comparable Rabi frequencies. This implies that the quantum state measured after 
propagation can be completely different from the initial state, even though the mean values of the 
field are unaltered. 

PACS numbers: 42.50.Gy,42.50.Lc,42.50.Ar,42.50.Pq 



Electromagnetically induced transparency (EIT) [ij is 
a technique that can be used to eliminate fluorescence 
from an atom illuminated with light whose frequency is 
equal to a particular atomic transition. This technique 
can be used in systems of three-level atoms in A configu- 
ration 0] , see Fig. ^ In this configuration a mode of the 
field, called the pump field, interacts resonantly with one 
dipole transition, while another mode, the probe field, in- 
teracting with the second dipole transition, is tested for 
transparency. The linear response of the absorption of 
the probe field by the medium is described by the imag- 
inary part of the electric susceptibility. In Fig. [3 we plot 
the susceptibility as a function of the probe frequency 
(solid line). The maximum absorption of the probe field 
by the medium depends on the Rabi frequencies associ- 
ated with each atomic optical transition. The maximum 
occurs for a detuning from resonance which increases 
monotonically with the Rabi frequencies. 

Many recent works have investigated if this trans- 
parency, originally studied for classical fields, preserves 
the initial quantum properties of the probe field. For 
a classical pump field with Rabi frequency much larger 
than that of the probe, Lukin et al. [3j showed that the 
medium is transparent for the quantum state. Further- 
more, they demonstrated a transfer of the initial quan- 
tum state from the probe field to the atoms and in a 
second stage back to the field by varying the Rabi fre- 
quency of the pump laser. They proposed using this tech- 
nique as a quantum memory device. When both fields 
are treated quantum mechanically, Dantan et al. stud- 
ied the noise spectrum of the quadratures, when only the 
coherent pump field drives the atoms and the probe field 
is initially in a broad-band squeezed vacuum. If the fre- 
quency equals that of the atomic transition, the medium 
is transparent. For other frequencies, there is absorp- 
tion of the quantum properties. The spectral absorption 
varies in a similar manner as the transparency curve for 
the mean value of the field. The transparency for the 
vacuum squeezed state was confirmed experimentally by 



FIG. I: The atoms have a A configuration with stable or 
metastable states |f) and |2) and common excited state |e). 
The transitions \i) ^ |e) with dipole coupling constants Qj un- 
derlie spontaneous decay of rates 7-, (j = 1,2), the linewidth 
of |e) is 7 = 71 -I- 72. 



Akamatsu et al. 0. 

In this Letter, we discuss the propagation of a quantum 
state in a medium of three level atoms showing EIT in the 
stationary regime. We focus on the case of a squeezed ini- 
tial state of the probe field and a coherent pump field. We 
treat both fields quantum mechanically and in contrast 
to previous work we do not put any constraint on the 
probe Rabi frequency. An analytical solution is given for 
this general case. Our main result is that the absorption 
of the initial squeezing is accompanied by an oscillatory 
interchange of the quantum properties between the pump 
and probe field while traveling through the medium. This 
oscillatory behavior is present for frequencies where the 
absorption of the mean value of the field is usually negli- 
gible. The effect is maximally pronounced for the case of 
equal Rabi frequencies. As a consequence, the probe and 
pump states after propagation can be completely different 
from the input state. 

To probe the behavior of quantum fluctuations of light 
passing through an EIT-medium, we use squeezed states 
defined by |q;;^) = D{a)S{£)\0) with the squeezing oper- 
ator 5(f) = exp[(fa^^ — f*a^)/2] and displacement oper- 
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for the field operator a(z, t) gives 
d 



d_ 

dt 
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FIG. 2: Absorption spectrum of the probe's mean value 
(solid line) and the function P{uj,0), Eq. IIUII . occurring in 
the spectrum of quantum fluctuations (dashed line), in ar- 
bitrary units. H[ere, cu — refers to a resonant driving of 
the probe transition. Both curves are normalized to have the 
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and tjJnax'^ ~ +fl2 are slightly shifted. (Parameters: 

^1 



1,91 = 52.) 



where we have introduced the coupling constant gj = 
pEj/h. Using the inversion operators vjj = a^e — crjj, 
the Heiscnberg-Langevin equations of the atomic degree 
of freedom for resonant driving read 

d 1 

— TOi ==-(-71 - 7)(1 + mi+ W2) - 2igi{(Teiai - ajcrie) 

- ig2{<Je2a2 ~ al(J2e) + /^i , 

d 1 

'gl^2 =^(-72 - 7)(1 + ^71+ K72) - igi{aeiai - a\aie) 

- 2ig2(CTe2a2 - 4^26 ) + /roa ' 

TrrCTie = - -o-ie + igiTUiai - ig2cri2a2 + he, 
ot 2 

^ _L • • ^ f 

^Cr2e = - 2'^2e + ig2^2a.2 ~ «5lO'2iai + J2e, 

d . t . 



ator D{a) 



exp[aa 



i], where a and at are the ^^'^21 - - maicr2e + «<?2CTeia2. 



(2b) 



annihilation and creation operator of the mode under 
consideration. After propagation along the z-axis, we 
analyze the fluctuation spectrum 



S{co) 



(1) 



in the steady state of the quadrature fluctuation SY^ (i) = 
6a{t) exp(— i0) + 6a'^{t) exp(i0), where 6a — a — (a). 

We consider two quasi monochromatic one dimen- 
sional beams of light propagating along the z axis in 
a medium of N three-level atoms in A-configuration. 
The excited state |e) of the atom with total linewidth 
7 = 71 -|- 72 can decay spontaneously with rate 7^ 
into the lower electronic state \j) (j = 1,2). To de- 
scribe the propagation of the two beams, we use a mul- 
timode representation of a pair of electromagnetic fields 
Ej = £jaj{z,t) exj>[iki^_jZ — uJhjt] + h.c. and treat the 



medium in a continuum approximation |J,|y{. Here, \£j\ 
is the vacuum electric field at the laser's carrier frequency 
WL,j = ckj^j, and aj{z,t) is the envelope operator of the 
corresponding field j, which is slowly varying in space 
and time. For the medium we assume the inter-atomic 
distance to be much smaller than the shortest relevant 
wave length of the laser light and introduce the continu- 
ous atomic operators iJi^^{z) = limAz-^o jvSI S 

is the individual atomic operator 
of atom j at position zj . 

The atoms interact with the electric fields via their 
dipole moment p; we treat this interaction in the rotating 
wave approximation. Then, in the slowly varying enve- 
lope approximation, the Heisenberg equation of motion 



0) 

fJ.D , 



where cr^i) — 



This system of equations (|2a|l and Ij2b|l is easily inter- 
preted: The polarization field aje{z,t) serves as a source 
for the electric fields, whereas the propagating light in 
turn drives the atomic media via the dipole interac- 
tion terms (x gj. The fj are delta-correlated, collec- 
tive Langevin operators which account for the noise in- 
troduced by the coupling of the atomic system to the 
free radiation field. They have vanishing mean values 
and correlation functions of the form {fx{z, t)fy{z, t')) = 
jfDxy5{t — t')5{z — z'). The diffusion coefhcients D^y can 
be obtained from the generalized Einstein equations 0. 
For this particular problem they are listed in ^ 

We will use the common technique 3 0| of 
transforming Eqs. (|2a|l and (|2b(l into stochastic c- 
number equations which serve to calculate correla- 
tion functions of the atomic and field operator up 
to second order. To fix the corresponding order 
of operators we use the "normal" order convention 
a^, a|, (Te2, cTei, 0-12, roi, 1372, (721, die, cr2e, ai, 02. All Calcu- 
lated results from the c-number equations have to be con- 
sidered to represent the results using operators in this 
order. The equations of motion for the c-number quan- 
tities are equivalent to those for the operators, except 
that a redefinition of the diffusion coefficients is neces- 
sary. These new diffusion coefficients are listed in 0. In 
the following, we use the same symbols for operators and 
their c-number substitutions, except that for the field op- 
erators Oj we instead use aj. 

In order to solve the c-number counterparts of 
Eqs. (|2a|l and (|2b|l . it is convenient to collect the sys- 
tem quantities and the fluctuation forces in vectors 
yi^{z,t) = (q;2, a*, cre2, . . . , cr2e, ai, 012) and f'^(z,t) = 
(0, 0, /e2, . . . , /2e, 0, 0), whcrc for the components we 
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choose the same order convention introduced above. For 
a large number of atoms, it is reasonable to assume a 
steady state for the mean values (x) with small fluc- 
tuations (5x, i.e. we write x = (x) -I- Sx. This allows 
us to treat the problem perturbatively for small fluctu- 
ations (5x oc 0(1/ VN) 11 [13. The steady state val- 
ues (x) cx 0(1) are found from the equations of motion 
()2a|l and l|2b|l to zeroth order in Sx, after setting all time 
derivatives to zero. Due to EIT the steady state values 
are z-independent. To first order, the system of equations 



^Sx{s, t) = M{s)5x{s, t) + f (s, t) + g(i) (3) 
at 

is linear in the fluctuations, where (5x(s, t) ~ 

00 

J exp[—sz]6x{z,t)dz (and analogously for f) denotes the 


Laplace transform of 5x{z,t). The elements of the 

matrix A4{s) are calculated from the equations H2a|) 

and Ij2b|l to first order in 6x. In g(t) we collect 

the initial conditions at z = following from the 

Laplace transform of Eqs. H2a(l . Equation (jSJ can be 

cast into algebraic form using the Fourier transform 
00 

Sx(s,u)) — J dt exp[~iujt]5x{s,uj) /\/27r. From the so- 

—00 

lution 6x{s,uj) = [A4 +iuj]^^{i + g) we can construct the 
correlation matrix 

{6x{s,uj)Sx\s',oj')) ^6{LU + iu')x 

[M{s) + liu]-^ {V{s, s') + g{iu)) [M^{s') - lu,]-^ 

(4) 

with VS{uj + uj') = (fft), gS{uj + uj') = (ggt) and (fgt) = 
(gf^) = 0. For fields entering the medium in a squeezed 
state with real squeezing parameter (j ~ 1,2), we have 



—c^Sij cosh^i sinh^i , 
c^Sij sinh^ , 



(5) 
(6) 



with all other coefficients of Q vanishing. 

We can now use the result Q to calculate the spectrum 
of fluctuations of the electric fields. To this end we recall 
the relation 

{SYf{z,Lu)SYf{z,uj')) - S{uj+lu')[S,{lu) - 1] , (7) 

following from the Wiener Khinchine theorem, connect- 
ing the fluctuation spectrum Sj{uj), Eq. with the 
correlation of the quadrature fluctuations SYj^{z,u)) = 
Saj{z,uj)e~'^^ + Sa*{z, —uj)e^^. Recall that the c-number 
result represents the normal ordered spectrum neglect- 
ing vacuum fluctuations. To take these into account, 
we write [Sjioj) — 1] in Eq. Q in order to ensure that 
Sj{uj) — 1 for a coherent state. The spectrum Sj{uj) 
can now be calculated by evaluating the left-hand side 
of Eq. Q with the help of the corresponding matrix 



elements from Eq. Q using a two-dimensional inverse 
Laplace transform in s and s' . 

As initial conditions we specify ^1 = and real ^2 = 
that is, the probe held (j = 2) is in a broad band 
squeezed state, whereas the pump field is coherent. All 
frequency components of both pump and probe beam, are 
in a squeezed/coherent vacuum {aj = for to ^ wl), only 
the carrier frequencies of the two beams are displaced by 
real a^, thus driving resonantly the atomic transitions 
with Rabi frequency Qj = \gjaj\. 

Using these initial conditions, we find for the fiuctua- 
tion spectra of pump {j — 1) and probe {j = 2) 



Si{z,uj) ==1 



/(e,g) 

^4 



(8) 



2n|rj2 



with the resonance curve 



P(c.,A) 



(7/2)2(^2 + (cj2_f]2)2 ■ 



(10) 

Furthermore, we have defined /(^, 6*) = 1 — e'^^ cos^ 9 — 
e"2f sin^ 9 and n'^ = nj + fl^. Indeed, for z = 0, we find 
iSi(0,a;) = 1 and 82(0, uj) — e^^ for the = quadrature. 
In Fig. 13 (a) we have plotted the fluctuation spectrum of 
the probe pulse as a function of frequency and propaga- 
tion length z. In parts (b) and (c) of Fig.Ocuts along the 
Lo and z direction are shown for both, pump and probe 
fields. 

We start the discussion of Eqs. © and © with ana- 
lyzing their asymptotic behavior. For z — > 00 we obtain 

02o2 

5i(z,c.)«i-/(e,0)^, 

52(z,c.)«l-/(e,0)|^. 

The distance where this asymptotic behavior is 
achieved is governed by the exponentials in Eqs. (|HJ| and 
0, and is of the order of Zabs = l/7-P(w,0). When 
ill 7^ 0, the asymptotic quadrature fluctuations of the 
probe show absorption of the initial squeezing of the 
field. As P(u;,0) follows qualitatively the probe field 
mean value transparency curve, see Fig. |21 so does the 
distance where this asymptotic behavior takes place. The 
pump field shows also squeezed fluctuations in this limit, 
but from the analysis of different quadratures 9 it fol- 
lows that the state is no longer a minimum uncertainty 
state. For equal Rabi frequencies Oi =02, both fields 
have asymptotically the same fluctuations. Similar cor- 
relations are known in cavity EIT and from the effect 
of pulse-matching 0, ^| . 
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FIG. 3: Quantum properties of laser beams propagating through an EIT medium, (a) The fluctuation spectrum of the probe 
beam, initially (at z = 0) in a squeezed state with squeezing parameter ^ = —3 propagates along the z-direction. Shown is 
the variance of the quadrature — for different frequencies ui/j and positions zC/7, where C is the prefactor of Eq. llUt . 
(b) The fluctuations of pump (dashed) and probe (solid) as a function of position for fixed frequency u = 0.257 and 6^ = 0. 
For the probe, the fine dotted line shows a simulation with finite lifetime 5OO/7 of the ground state coherence, c) Fluctuation 
spectrum of pump and probe at fixed position z = IOO7/C and G = 0. (Parameters: — = 7, 5i = 52 = 7/6O, — —3.) 



In the case of strong asymmetric driving Q2 ^ ^^i, it 
follows from Eqs. ||SJ) and © that the pump field always 
stays in a coherent state while the initial fluctuations of 
the probe field are exponentially damped during propa- 
gation through the medium. The intensity of the absorp- 
tion is given by ^P{uj, 0). As already mentioned, P{uj, 0) 
behaves like the mean value transparency curve, so that 
we find a comparable behavior of the absorption of the 
probe's quadrature fluctuations as for its mean value ab- 
sorption. Thus, in this limit, our result reproduces the 
analysis reported in Ref. 4]. 

When f2i « i72, however we flnd an entirely novel be- 
havior. In order to bring out this new aspect most clearly, 
we assume 7-P(w, 0)z ^ 1, that is, we consider positions z 
where the exponential absorption of the quadrature fluc- 
tuations can be neglected. Then, Eqs. © and can 
be approximated as 

S,{C,u)^l-f{^,e)^^B\n'C (11) 

52(C, -) - 1 - { {^l mini cos^ c} 

(12) 

with ^ = z/zqsc, where we introduce the oscillatory 
length scale Zqsc — 2/ P{LU,n)uj. Equations Hll() and 
(|12|l clearly display the oscillatory transfer of the ini- 
tial quantum properties of the probe to the pump and 
back while traveling through the medium. During this 
oscillatory behavior, the pump and probe field stay in a 
minimum uncertainty state. Also the sum of the fluctu- 



ations Si + S2 is conserved. When fli = fl2, this oscilla- 
tory transfer is maximal and all the fluctuation proper- 
ties oscillate between the probe and the pump field dur- 
ing propagation in the medium. The length scale of the 
oscillatory transfer Zqsc can be much smaller than the 
absorption length scale Zabs- In Fig- El the interplay of 
both scales, oscillatory and absorption, can be clearly 
observed. 

The oscillatory behavior implies that the outgoing 
quantum field can be completely different from the in- 
coming field, although the mean values stay exactly the 
same. In other words, the medium is not transparent for 
the quantum properties of the field, at least not in the 
sense that the quantum state can traverse the medium 
unaltered. The oscillatory behavior is qualitatively dif- 
ferent from the absorption and is "coherent" since it does 
not imply loss of quantum properties. 

To check if these oscillations still persist with a damp- 
ing rate of the ground state coherence we performed a 
numerical simulations. For a typical value — I/5OO7 
the resulting curve (dotted) shown in Fig. O b) almost 
coincide with the ideal case. A complete discussion of 
the effect of decoherence will be published elsewhere. 

Present technology allows for an experimental obser- 
vation of this effect. The initial state can be constructed 
by mixing a laser tuned to the probe optical transition 
with a wide band squeezed field which covers the trans- 
parency window. With a Rabi frequency Qi = CI2 = "f 
and an observation frequency to « O.I7, Zabs is ten times 
Zosc- For a density p « 10^, we have a maximum transfer 



of quantum properties between the fields for z « 7cm, 
while Zabs ~ 70cm. When w ~ 0.257 we have a maxi- 
mum transfer of quantum properties between the fields 
for z K, 4.5cm {zC ^ 25 in Fig. Ob)) while ^abs ~ 12cm 
{zC /7 « 64) . Parameters on this order are found in var- 
ious BIT experiments [sLIllllst. 

In conclusion, we have shown that in the propagation 
of pump and probe beams through an EIT medium, the 
quantum states are not conserved, except for the carrier 
frequencies which drive the atoms on two-photon reso- 
nance. Apart from absorption of quantum fluctuations, 
which takes place in approximately the same frequency 
range as the absorption of the mean values, we found 
a novel characteristic behavior, which is most strongly 
pronounced if the two beams have comparable Rabi fre- 
quencies, consisting of an oscillatory transfer of the ini- 
tial quantum properties between the probe and the pump 
field. This effect could be observed in current state-of- 
the-art experiments. 
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